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The characteristic function of soHton phase jitter is found analytically when the 
soliton is perturbed by amplifier noise. In additional to that from amplitude 
jitter, the nonlinear phase noise due to frequency and timing jitter is also 
analyzed. Because the nonlinear phase noise is not Gaussian distributed, the 
overall phase jitter is also non-Gaussian. For a fixed mean nonlinear phase 
shift, the contribution of nonlinear phase noise from frequency and timing 
jitter decreases with distance and signal-to-noise ratio. 
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1. INTRODUCTION 

The phase jitter of soliton due to amplifier noise, like 
Gordon-Haus timing jitteri, is usually assumed to be 
Gaussian distribute d^i'^i'^i^ . When the phase jitter of 
soliton was studied, the phase jitter variance was given 
or measured and the statistics of soliton phase is not 
discusse d^i?'^'Ti? . 

For non-soliton systems, the statistics of nonlinear 
phase noise is found to be non-Gaussian distributed 
both experimentallji^ and theoreticalljiiSiiiiiSii^. How- 
ever, those studieaSiiSiiiiiSii^ just includes the Gordon- 
MoUenauer effeclii^ that is the nonlinear phase noise 
induced by the conversion of amplitude to phase jit- 
ter due to fiber Kerr effect, mostly self-phase modula- 
tion. Based on the well-developed perturbation theory 
of soliton^ii^iiSiii, phase jitter can also be induced by 
the interaction of frequency and timing jitter. In this 
paper, the statistics of the soliton phase is derived in- 
cluding the contribution of timing and frequency jitter 
induced nonlinear phase noise. The characteristic func- 
tion of soliton phase jitter is derived analytically, to our 
knowledge, the first time. The probability density func- 
tion (p.d.f.) is simply the inverse Fourier transform of 
the corresponding characteristic function. 

Most optical communication systems use the 
intensity of the optical signal to transmit infor- 
mation. Direct-detection differential phase-shift 
keying (DPSK) signaling has renewed attention 
recently^^i^^i'^°i^^i^^i^'^i'^'^i^^i^^i^^, mostly using return- 
to-zero (RZ) pulses for long-haul transmission and 
encode information in the phase difference between 
two consecutive pulses. To certain extend, a soliton 
DPSK system may be a good approximation to phase 
modulated dispersion managed solitonA or RZ signal. 
With well-developed perturbation theorjiSii^*iSiii, the 



distribution of the soliton phase jitter can be derived 
analytically. 

The error probability of DPSK soliton signal was cal- 
culated in Ref. |0 using the method of Refs. andlsol 
without taking into account the effect of phase jitter. If 
the phase jitter is Gaussian distributed, the system can 
be analyzed by the formulas of Ref. |^ The phase jitter 
may be indeed Gaussian distributed in certain regimes 
around the center of the distribution^-^^, especially if the 
p.d.f. is plotted in linear scale. The tail probability less 
than, for example, 10^^, is certainly not Gaussian dis- 
tributed. As optical communication systems are aimed 
for very low error probability, a careful study of the statis- 
tics of the soliton phase is necessary to characterize the 
performance of the system. 

The remaining parts of this paper are organized as 
following: Sec. |21 gives the stochastic equations of the 
phase jitter according to the first-order soliton perturba- 
tion theory; Sec. 13 derives the characteristic function of 
soliton phase jitter; Sec. 0]presents the numerical results; 
and Sees. [51 and are the discussion and conclusion of 
the paper, respectively. 

2. STOCHASTIC EQUATIONS FROM SOLI- 
TON PERTURBATION 

From the first-order perturbation theory, with amplifier 
noise, the soliton parameters evolve according to the fol- 
lowing equations^ii^iiSiii 
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where 5f{ } and Q{ } denote the real and imaginary part 
of a complex number, respectively, n((^, r) is the amplifier 
noise with the correlation of 



i?{n(Ci,Ti)n(C2,r2)} =a25(Ci -C2)<5(ti -T2), (5) 

A{C,), ri(C), T{C,), and are the amplitude, frequency, 
timing, and phase parameters of the perturbed soliton of 



go(r,C) = A(C)sech{A(C)[T-T(C)]} 

X exp [-iniOr + i(^(C)] (6) 

with initial values of A{0) = A and Q{0) = 0(0) = 
r(0) = 0. Functions related to soliton parameters are 
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Using Eqs. ©, lfTT|) . and lfTC|) . the phase jitter is 



(18) 




wr2(C2)dC2 +u't(Ci)J dwn(Ci) 

(19) 



where is a zero-mean Wiener process with autocorre- 
lation function of 



fn = tanh[A{T-T)]q*, 

JT = —4—90. 
1 



(7) 

(8) 

(9) 



{1 _ A{t - r)tanh[A(T - T)]} q*. (10) 



From both Eqs. ^ and (0), we get 

A{C) = A + WAiO 
^{0 = wniO 



(11) 
(12) 



where wa and wq are two independent zero-mean Wiener 
process with autocorrelation functions of 



E{wa{Ci)'wa{C2)} 
E{wniCi)wniC2)} 



cr^ min(Ci,C2), (13) 
4min(Ci,C2), (14) 



where a\ — Aa^ and ct^ = Aa^/S^iii. Defined for the 
amplitude, the signal-to-noise ratio (SNR) as a function 
of distance is 



A^ 



A 



Using Eqs. Q and ifT^ . the timing jitter is 



(15) 



T(C) 



(16) 



where wt is a zero-mean Wiener process with autocorre- 
lation function of 



E{wT iCi)wTiC2)} = o-rmin(Ci,C2) 



(17) 



with 



5.7.17 



E{w4,{Ci)w4,iC2)} = cr0niin(Ci,C2) 



witl 
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(20) 



(21) 



The Wiener processes of wa,wq,wt, and are in- 
dependent of each other. The amplitude [Eq. Hll() ]. fre- 
quency [Eq. ifT^ ]. and timing [Eq. ((TE|) ] jitters are all 
Gaussian distributed. From Eq. H19|l . it is obvious that 
the phase jitter is not Gaussian distributed. If Eq. Q 
is linearized or all higher-order terms of Eq. H19|l are ig- 
nored, the phase jitter is Gaussian distributed and equals 
to 0(C) « A J^WA{Ci)dCi + w^{() [Ref.Q. The charac- 
teristic function of the phase jitter Eq. H19|l will be de- 
rived later in this paper and compared with Gaussian 
approximation. 



3. CHARACTERISTIC 
PHASE JITTER 



FUNCTIONS OF 



In the phase jitter of Eq. (|19|) . there are three indepen- 
dent contributions from amplitude jitter (the first term), 
frequency and timing jitter (the second and third terms), 
and the projection of amplifier noise to phase jitter w^. 
In this section, the characteristic functions of each indi- 
vidual component are derived and the overall character- 
istic function of phase jitter is the product of the charac- 
teristic functions of each independent contribution. 

A. Gordon- Mollenauer Effect 

The first term of Eq. (|19|1 is the Gordon-MoUenauer 
effecOi of 



0gm(C) 



[A + WA{(:i)?dCi 



(22) 



induced by the interaction of fiber Kerr effect and ampli- 
fier noise, affecting phase-modulated non-return-to-zero 
(NRZ) and RZ signaU^iiJ^. 
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The characteristic function of Gordon-Mollenauer non- 
linear phase noise is given bjiii^i^ 



*0gm(C)('') = sec* (^C(^aVp) 
2aA 



X exp 



A" I— ( ^ 

- — VJ^tan \CaA 



(23) 



The above characteristic function Eq. 1)23(1 can also be 
derived from Eq. \K.1\ of the appendix. 

The mean and variance of the phase jitter Eq. (|22|l are 



<0Gm(C)>=-J ^*0gm(C)('^) 



v=0 



" ^2 /■2 



(24) 



and 



^«n(Ci)dCi 



u;T(Ci)dwn(Ci 



-wn(C) / «^o(Ci)dCi- 
Jo 



(27) 



From Eq. ((A. 12(1 of the appendix, the characteristic func- 
tion of 0a,T(C) is 



(28) 



The mean and variance of the phase jitter of Eq. 1(26(1 



are 



-^4c^ (29) 



and 



>gm(C) 



^^0gm(C)('^) 



<0gm(C)>' 



1 
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(25) 



respectively. The first term of Eq. 1(251) increases with C'^, 
conforming to that of Ref. El Given a large fixed SNR 
of A^/{a\C,) [Eq. (O], the second term of Eq. ^ is 
much smaller than the first term and also increases with 
C^. Note that the first term of the mean of Eq. ((24(1 is 
also larger than the second term for large SNR. 

The characteristic function of Eq. ((23(1 depends on two 
parameters: the mean nonlinear phase shift of j4-^C/2 and 
the SNR of Eq. ((15(1 . Given a fixed mean nonlinear phase 
shift of A^C/2, the shape of the distribution depends only 
on the SNRii. 

Based on Eq. p9(l . comparing Eq. 1(23(1 with the non- 
soliton case of Ref. .11. , the mean and standard devia- 
tion of the Gordon-Mollenauer phase noise of soliton are 
about half of that of non-soliton case with the same am- 
plitude A as the NRZ or RZ leveiii. 

B. Frequency and Timing Effect 

The frequency and timing jitter contributes to phase jit- 
ter by 



M'o(C2)dC2du)n(Ci) 



JO 

C 



u;T(Ci)dwn(Ci) 



(26) 



as the second and third terms of Eq. II19() . 

By changing the order of integration for the second 
term of Eq. ((25(1 . we get 



>n,T(C) 



= - ^*0o,t(C)('^) 



- <0f2,T(C)>^ 



(30) 



respectively. 

Comparing the means of Eqs. ((24(1 and II29() . in terms 
of absolute value, the mean nonlinear phase shift due to 
Gordon-Mollenauer effect is much larger than that due to 
frequency and timing effect. Comparing the variances of 
Eqs. ((25(1 and l(3U() . the variance of nonlinear phase noise 
due to Gordon-Mollenauer effect is also much larger than 
that due to frequency and timing effect. 

Unlike the Gordon-Mollenauer effect, the characteris- 
tic function of Eq. 1(28(1 . from the appendix, is not deter- 
mined only on the SNR and the mean nonlinear phase 
shift Eq. (E3l. 

C. Linear Phase Noise 

The last term of Eq. ((19(1 gives the linear phase noise of 



0ln(C) = 
with a characteristic function of 



*0ln(C)(^) = exp 



1 2 /• 2 



(31) 



(32) 



From the characteristic function of Eq. I(32() . the linear 
phase noise depends solely on the SNR [Eq. ((T3|) ]. 

The characteristic function of the overall phase jitter 
is the multiplication of the characteristic functions 
of Eqs. (ESJ, (EHl), and 

Although the actual mean nonlinear phase shift is 

<0(C)> = <0f2,T(C)> + <0GM(C)>, (33) 

we mostly call A^C/2 the mean nonlinear phase shift as 
a good approximation in high SNR. 
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4. NUMERICAL RESULTS 

The p.d.f. is the inverse Fourier transform of the cor- 
responding characteristic function. Figures ^ show the 
evolution of the distribution of the phase jitter [Eq. H19(l ] 
with distance. The system parameters are A = 1 and 
cr^ = 0.05. Those parameters are chosen for typical dis- 
tribution of the phase jitter. 

Figures n^a), (b), (c) are the distribution of Gordon- 
MoUcnaucr nonlinear phase noise [Eq. (|23|l ]. frequency 
and timing nonlinear phase noise [Eq. (|28|) ]. and the lin- 
ear phase noise [Eq. (|32() ]. respectively, as components of 
the overall phase jitter of Eq. H19|l . Figure^^d) is the dis- 
tribution of the overall phase jitter Eq. fll^ . The p.d.f.'s 
in Figs, n are normalized to a unity peak value for illus- 
tration purpose. The x-axis of individual figure of Figs.^ 
does not have the same scale. From Figs. Q] the nonlin- 
ear phase noises from Gordon-MoUenauer effect and fre- 
quency and timing effect are obvious not Gaussian dis- 
tributed. With small mean and variance, the nonlinear 
phase noise from frequency and timing effect has a very 
long tail. 

Figures|3plot the p.d.f.'s of Figs.^in logarithmic scale 
for the cases of C = 1,2. The Gaussian approximation is 
also plotted in Figs. 12 for the overall phase jitter 
In both cases of ^ = 1,2, the Gaussian approximation is 
not close to the exact p.d.f.'s in the tails. However, if the 
p.d.f.'s are plotted in linear scale, Gaussian approxima- 
tion may be very close to the actual distribution, espe- 
cially for large phase jitter—. The p.d.f's in Figs. |21 are 
not normalized to a unity peak. 

From both Figs. Q and [3 the nonlinear phase noises 
of 0GM and (f>n^T are not symmetrical with respect to 
their corresponding means. While (f>QM spreads further 
to positive phase, 4'n,T spreads further to negative phase. 
Plotted in the same scale, the nonlinear phase noise of 
(f>GM due to Gordon-Mollenauer effect is much larger than 
the nonlinear phase noise of 4>n,T due to frequency and 
timing effect. 

The p.d.f.'s in Figs. ^ cannot cover all possible cases. 
While both the Gordon-Mollenauer and linear phase 
noises depend on the mean nonlinear phase shift A^C^jl 
and SNR, the nonlinear phase noise induced by frequency 
and timing effect does not have a simple scaled relation- 
ship. 

For a mean nonlinear phase shift of \^C, — 1 rad^*. 
Figures 01 plot the distribution of the overall phase jitter 
[Eq. (dSJ] for a SNR of 10 and 20 for C = 1, 10. After a 
scale factor, the distributions of both Gordon-Mollenauer 
and linear phase noise are the same as that in Figs. H 
In additional to the overall phase jitter. Figures |2I also 
plot the distribution of the nonlinear phase noise from 
frequency and timing effect of 4)a,T- 

For a fixed mean nonlinear phase shift and SNR, from 
Figs. 13 the nonlinear phase noise from frequency and 
timing effect of (pn^riO has less effect to the overall phase 
jitter for long distance than short distance. Figures^are 
plotted for short distance of C < 3 to show the contri- 
bution of frequency and timing jitter to nonlinear phase 



noise. The effect of <?!>n,T(C) is smaller for large SNR of 
20 than small SNR of 10. The main contribution to the 
overall phase jitter is always the Gordon-Mollenauer ef- 
fect and the linear phase noise. 

5. DISCUSSION 

The phase jitter of Eq. (|19|l is derived based on the first- 
order perturbation theorySii5iiS*ii of Eqs. (P) to (@J. The 
non-Gaussian distribution is induced by the higher-order 
terms of Eq. H19() or the nonlinear terms of Eq. Q). 
Second- and higher-order soliton perturbation'^^'^^ may 
give further non-Gaussian characteristic to the phase jit- 
ter. Currently, there is no comparison between contri- 
butions of the higher-order terms of Eq. Q and higher- 
order soliton perturbation. 

In this paper, like almost all other 
literaturesiiSi^iiSiiiiiSiiSiii, the impact of amplitude 
jitter to the noise variances of cr^, cr^, a^, and cr^ is 
ignored. The noise variances of CTq, a^, and cr^ 
are assumed independent of distance. If the amplitude 
noise variance is a\ = A{()it'^ with dependence on the 
instantaneous amplitude jitter, amplitude, frequency, 
and timing jitters are all non-Gaussian distributed"^^. As 
an example, amplitude jitter is non-central chi-square 
distributedi2^i2&. However, the statistics of phase jitter 
[Eq. (|19l) ] does not have a simple analytical solution 
when the noise variance depends on amplitude jitter. 
With a high SNR, the amplitude jitter is always much 
smaller than the amplitude A{0) = A. Even in high 
SNR, the phase jitter is non-Gaussian based on Eq. H19|l . 

6. CONCLUSION 

Based on the first-order soliton perturbation theory, the 
distribution of soliton phase jitter due to amplifier noise 
is derived analytically the first time. In additional to the 
main contribution of Gordon-Mollenauer effect, the non- 
linear phase noise due to frequency and timing jitter is 
also considered. Induced by Gordon-Mollenauer effect or 
frequency and timing jitter, the nonlinear phase noises 
are not Gaussian distributed, neither does the overall 
phase jitter. For a fixed mean nonlinear phase shift, the 
contribution of nonlinear phase noise from frequency and 
timing jitter decreases with distance and SNR. 

Appendix A 

Here, we find the joint characteristic function of 

VI = f^wUCi)dCu (A.l) 
Jo 

V2 = f WT{(:i)dwn{Ci), (A.2) 
Jo 

V3 = wniO I wn{Ci)dCi. (A.3) 
Jo 
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(a) 0gm(C) 



(b) <^o,t(C) 



3 -1.5 

(c) <^ln(C) 




3 -1 

(d) <^(C) 

Fig. 1. The distributions of soliton phase jitter for difference distance for A = 1, af^ ~ 0.05. The distributions are 
normahzed for a unity peak. The x-axis is not in the same scale. 

By changing the integration order, we get 



10 Jo Jo 



(A.4) 



The joint characteristic function of (pi, (p2, and (ps is 



*Vl,V2,y3(^l'f'2,l'3) = £^{exp(jj/i(/3i +jV2(/32 +ji^3V3)}- (A. 5) 

Similar to option pricing with stochastic volatilitji^l, the expectation of Eq. (jA.Sp can be evaluated in two steps, first 
over wt and than wq. In the average over wt, it is obvious that ip2 is a zero- mean Gaussian random variable with a 
variance of cr| [wn{C) - W(i{Ci)f dCi, we get 



t2 ,,2 /-C 



K7(C)-^o(Ci)] dCi 



+3^1 i Wo(Ci)dCi + jV3Wn(C) / wo(Ci)dC; 



wm) + (jV3 + 4^^2)wo(C) / wniQMQ 



t2 7/2 \ /-C 



«^o(Ci)dCi 



(A.6) 
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First of all, we liavo^9i^?i??i?? 




where uji^2 = {^1,^2)'^ and 
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(To tan ( VJiJ^anC) 



1 

jW3 



sec 



J ^3 
1 

i^3 



sec (v/ji^o-nC) - 1 
tan (Vj^^ctoC) 



As a verification, if approaches zero, the covariance matrix is 



lim Cljuj^) 



c 

2*5 3'' 



that is the covariance matrix of the vector of 



(A., 



(A.9) 



(A.IO) 



withou t any dependence on the random variable (pi. Note that the equation corresponding to Eq. ljA.7p in Refs. 
|l(f39l does not have the hmit of E q. (|IT9)l . 

The characteristic function of Eq. (|A.7p is that of a correlated two-dimensional Gaussian random variable of with 
dependence to (pi. The first two terms of Eq. (|A.6|I is a quadratic (or bilinear) function of w^;, i.e., ^w^Al(jV2, jV3)w^, 
where 



J 1^3 +a^V2 



(A.ll) 



The characteristic function of the quadratic function of zero- mean Gaussian random variables is det[I — CM] 2 
[Ref.lll, where det[ ] denotes the determinant of a matrix. 
The joint characteristic function is 



sec 2 i^2jvi - a'^vlanC, 



det [I-C {2jvi - a^vl)M{jv2,iv3)Y- 
I 



(A.12) 



where T is the identity matrix. The substitute of joj^ 
by 2jvi — a\v\ is obvious by comparing Eqs. (|A.6|I and 

We can ge1>iii2£ 
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